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The heavy mesons consisting of a light and a heavy 
quark form an interesting laboratory for the study of 
QCD, the theory of the strong interaction. They are 
a close analogue to the hydrogen atom and play a similar 
role for the study of the QCD as hydrogen for quantum 
electrodynamics. According to heavy quark effective the- 
ory (HQET) fH, in the limit of infinite mass of the heavy 
quark, the heavy quark decouples from the degrees of 
freedom of the light quark. For orbitally excited states 
{L = 1), the total angular momentum of the light quark 
is jq — 1/2 or jq = 3/2. Combining jq with the spin 
of the heavy quark, four states forming two jq doublets 
are expected. For an infinite mass of heavy quark, the 
four states are degenerate in mass. However, corrections 
due to finite mass lead to fine structure splitting between 
the two doublets and hyperfine structure splitting within 
each of the doublets. Following the standard scheme 
the states with jq = 1/2 are named B*q and Bsi, and 
the states with jq = 3/2 are named Bgi and B*2. Often 
these four states are referred to as B** . 

If kinematically allowed, all four B** states are ex- 
pected to decay dominantly to B*K, BK or both. The 
states of the doublet with jq — 1/2 decay through an S- 
wave transition and are therefore expected to have broad 
mass distributions. The states with jq — 3/2 decay 
through a D-wave transition and therefore are expected 
to have narrow mass distributions. In the following we 
focus on the narrow doublet. While the Bgi decays only 
into B*K due to conservation of spin and parity, the B*2 
can decay to B*K and BK. If the B*2 mass is near the 
B*K threshold, the decay to B*K will be strongly sup- 
pressed compared to the decay to BK due to the available 
phase space. 

Several theoretical predictions for the basic properties 
of the B** states are available ^3]. The predictions for 
the Bsi mass range from 5805 to 5891 MeV/c^, and for 
B*2 5820 to 5903 MeV/c^ with a mass difference between 
both states in the range 12 to 20 MeV/c^. The natural 
widths of the two states are expected to be of order 1 
MeV/c^ with strong variation with predicted mass. 

While there is already considerable information about 
D** mesons, the analogous particles in the charm sector 

0. 01: experimental knowledge about the B** mesons is 
minimal. First evidence for at least one of the B** states 
was found by the OPAL experiment Evidence for a 
single state interpreted as B*2 was seen by the Delphi 
Collaboration 6] and a preliminary observation of this 
state was reported recently by the D0 Collaboration 0] . 

In this Letter, we report on the observation of two 
states consistent with the jq = 3/2 doublet of the B** 
decaying to B+K- and B*+K~ with B*+ B+j, 
where the photon is not detected. Because of the miss- 
ing photon, the observed B^i peak is shifted downward in 
mass by the B*-B mass splitting of 45.78 ± 0.35 MeV/c^ 

1. B+ mesons are reconstructed in two decay channels. 



B+ J/^jK+ with J/iIj n+n- and B+ D \+ 
with D K+TT-. The use of a specific particle state 
implies the use of the charge-conjugate state as well. We 
use data collected by the CDF II Detector at the Fermi- 
lab Tevatron between February 2002 and February 2006 
corresponding to a total integrated luminosity of 1 fb~^. 

The components of the CDF II detector Q used for 
this analysis are the magnetic spectrometer and the muon 
detectors. The tracking system is composed of a silicon 
microstrip detector Q surrounded by an open-cell drift 
chamber (COT) flo'|. Both components are located inside 
a 1.4 T axial magnetic field. Muons are detected in planes 
of multi-wire drift chambers and scintillators [11] in the 
pseudorapidity range |ry| < 1.0, where 77 = — lntan(^^/2) 
and 6 is the polar angle measured from the proton beam 
direction. Hadron identification is crucial for distinguish- 
ing kaons originating from B** decays from other parti- 
cles. It is provided by a combination of the ionization 
energy loss in the COT and a measurement by a time-of- 
fiight system [l^ . 

A three-level trigger system is used for the online event 
selection. The level 1 trigger system includes the eX- 
tremely Fast Tracker (XFT) track processor fl3| which 
finds charged-particle tracks in the COT and measures 
their azimuthal angle around the beam direction and 
transverse momenta. In level 2, the silicon vertex trigger 
fl^ adds hits from the silicon detector to tracks found by 
the XFT to provide measurements of impact parameter. 
The level 3 system confirms the selections using a version 
of the offline event reconstruction optimized for speed. 

The dimuon trigger Q requires two tracks of opposite 
charge matched to track segments in the muon cham- 
bers, where the mass of the pair is consistent with the 
J/ip mass. The displaced-vertex trigger [3] requires two 
tracks with large impact parameters. Additionally, the 
intersection of the tracks has to be displaced from the 
interaction point and a minimum transverse momentum, 
the momentum component perpendicular to the proton 
beam direction, is required for each track. 

In both samples, we reconstruct B** candidates by 
combining B^ candidates with K'^ candidates. In 
the dimuon (displaced-vertex) trigger sample, we form 
J/tp — > M^/^^ ^ K^TT^) candidates and combine 
each J/ip [D ) candidate with a track assumed to be 
a kaon (pion), constraining the tracks to an appropri- 
ate decay topology to form a B^ candidate. At this 
stage hadron identification is not used. In order to im- 
prove the mass resolution, we consider the quantity Q de- 
fined as m{B+K-) - m{B+) - M^- , where m{B+K~) 
and m{B^) are the reconstructed invariant masses of the 
5+ K~ pair and the B^ candidate, and M^- is the 
known kaon mass [3]. The predicted B^i {B*2) state 
mass translates to the region < Q < 73 MeV/c^ 
(48 < Q < 131 MeV/c^). 

For the selection of candidates, we use a chain of two 
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neural networks based on the NeuroBayes package 
in each of the decay channels. In a first step, a neu- 
ral network in each channel combines topological, kine- 
matic, and particle identification quantities for the 
and its daughters to form a single discriminant between 
_B+ mesons and background. The most important quan- 
tities are the impact parameter of the B"*", the projection 
of the displacement of its reconstructed decay point from 
the beamline on the direction of its transverse momen- 
tum, the transverse momentum of the decay's pion 
(kaon), and its impact parameter. The neural networks 
are trained on two classes of events corresponding to the 
signal and background samples. In the B^ J/^pK^ 
channel, we use a pythia [17| simulation for the signal 
sample and experimental data from the B^ mass side- 
bands 5190 - 5240 MeV/c2 and 5320 - 5395 MeV/c^ for 

the background sample. In the B~^ 'D'\^ channel we 
use only experimental data to train the B^ neural net- 
work. We use candidates from a signal region between 
5240 and 5310 MeV/c^ in the invariant mass as signal 
sample and data from a mass sideband between 5325 
and 5370 MeV/c^ as background sample. The events 
from the B^ mass sidebands are used also as signal with 
negative weight to account for the background in the sig- 
nal region. Based on the neural networks, we select ap- 
proximately 31000 B'^ signal events in the J/tp decay 

channel and 27200 in the T^tt"*" channel. 

In a second step, we select B** candidates based on 
the number of candidates per event and on an additional 
neural network for each B~^ decay channel. These neural 
networks use the same inputs as used by the neural net- 
works to select B~^ mesons as well as their discriminant, 
and kinematic and particle identification quantities for 
the kaon track of the B** decay. The particle identifi- 
cation of the kaon is the most important variable, fol- 
lowed by the neural network discriminant of the B'^ and 
the pseudorapidity of the kaon. The number of candi- 
dates per event is not used in the neural network due to 
the difficulty of modeling fragmentation and hadroniza- 
tion in the production of heavy quarks. We select only 
those events with fewer than four candidates because a 
lower number of candidates provides a better signal-to- 
background ratio. This cut is fixed without looking to 
the experimental data, based only upon the above as- 
sumption. The B** neural networks are trained on a 
combination of simulated events, containing only signal, 
and data events in the Q range to 200 MeV/c^ for back- 
ground sample. The number of real B** mesons in the 
background-training sample is too small to affect signif- 
icantly the learning process of the neural network. In 
order to avoid possible mass biases, the simulated sig- 
nal events have the same B~^K~ mass distribution as 
the events used for background in the neural network 
training. The value of the cut on the neural network 
discriminant for the final selection is chosen to optimize 



NucNNdata, where Ndata (Nmc) is the number of the 
selected candidates in data (simulation) in the Q range 
60 to 70 MeV/c^. This range has been chosen based on 
the mass the of previously seen B*2 state and therefore is 
not biased with respect to the unobserved state. We ver- 
ify that the observed B** masses do not depend on the Q 
range used for cut optimization. The Q distributions of 
the selected candidates are shown in Figs.[lja) and[IJb) 
for the two trigger samples separately, and in Fig. [Ijc) 
added together. Two peaks are visible, centered near 67 
MeV/c^ and near 10 MeV/c^. The wrong-sign combina- 
tions (filled area in Fig. [1]) do not show any significant 
structure. 




°0 50 100 150 2C 

m(B^K>m(B^)-M [MeV/c^] 



FIG. 1: Distribution of Q = m{B+K-) - m{B+) - M^- for 
the B*" candidates with (a) B+ -* J/i^K+, (b) B+ 
and (c) both B^ channels combined. The dotted line shows 
the result of a fit with the sum of a background function and 
two Gaussians. The filled area shows the Q distribution for 
the wrong-sign combination B^K'^ . 
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The two peaks in data can be interpreted as the two 
jq = 3/2 states of orbitally excited Bg mesons. The natu- 
ral interpretation is that the peak near 67 MeV/c^ stems 
from the B*2 — > B^K~ decay while the peak near 10 
MeV/c^ stems from the decay Bgi — > B*^K~ . Revers- 
ing the assignment of the two peaks would result in a 
larger mass difference between BI2 and Bsi with Bsi 
being heavier, which would be opposite to other heavy 
quark mesons. 

To extract the mean Q values for the two peaks, we use 
an unbinned maximum likelihood fit. Each of the peaks is 
described by a Gaussian shape. We use a phenomenolog- 
ical function to describe the background without distin- 
guishing different types of backgrounds. The functional 
form of the background shape is — Qyf exp (^7Q), 
where [3 and 7 are free parameters. The fit has three 
free parameters for each of the Gaussians and two free 
parameters for the background. The fit to each data 
sample separately gives results consistent between the 
two B^ decay channels. Therefore, we combine the 5+ 
channels to perform the final fit. The projection of the 
fit on the full sample is shown in Fig. [IJc). From the fit 
we extract Q{Bsi) = 10.73 ±0.21 MeV/c^ and QfyBl^) 
66.96±0.39 MeV/c^^ with yields N{Bsi) = 36±9 events 
and N{B*2) = 95 ±23 events, where all uncertainties are 
statistical. The widths of the Gaussians are consistent 
with expected detector resolutions. 

Systematic uncertainties on the Q value may arise from 
the tracking and fitting procedures. The sources of un- 
certainty from the tracking are the uncertainty on the 
track error matrix, which enters through vertex fits, and 
the uncertainty on the material and magnetic field distri- 
bution inside the tracking volume. Based on a detailed 
study performed for the measurement of mass and width 
of the orbitally excited D** states [3], we assign a com- 
bined systematic uncertainty of 0.14 MeV/c^ due to the 
tracking effects. We study effects on the fitting procedure 
of the unknown background shape and the simplification 
by the single-Gaussian signal description. In both cases 
we generate a large number of samples of the same size as 
the data. For the background shape study we use a prob- 
ability density function proportional to a fit on data with 
a third order polynomial as background function instead 
of the default one for sample generation. In the signal 
shape study the sum of two Gaussians, which have the 
width one gets by fitting each decay channel separately, is 
used for sample generation. Each of the samples is then 
fitted using the default fit model and the pull distribu- 
tions are examined. In both cases the pull distributions 
are consistent with a Gaussian with a mean of zero and 
unit width. Therefore, we do not assign any systematic 
uncertainty arising from the fitting procedure. The re- 
sulting Q values are 

Q{Bsi) = 10.73 ± 0.21 (stat) ± 0.14 (sys) MeV/c^, 
Q{Bl2) = 66.96 ± 0.39 (stat) ± 0.14 (sys) MeV/c^ 



By adding the known values [3] of Mb* and Mj^- to 
Q{Bsi) and Mb+ with Mf^- to Q{B*2), we obtain 
m{Bsi) = 5829.4±0.7MeV/c2 andm(S;2) = 5839.6±0.7 
MeV/c^. The statistical and systematic uncertainties 
on the Q value and the uncertainties on the masses 
of and 5+ or B*^ are added in quadrature. Fi- 
nally, the mass difference of the two narrow i?** states 
is Ato(B*2,Bsi) = 10.5 ± 0.6 MeV/c^ where we add 
Mb- - Mb = 45.78 ± 0.35 MeV/c^ to the difference of 
the two measured Q values. 

To estimate the statistical significance of each of the 
two peaks, we repeat the fit without the term for one of 
the Gaussian peaks and again without the other peak. 
For each peak we form A£ = — 21n£o/'C, where C is 
the value of the likelihood function of the original fit and 
is the value for the fit without one of the peaks, and 
measure A£ = 48.7 (74.5) for B^i {B*2)- We generate 
samples with background according to the background 
function from the fit to the data and one of the peaks. 
For each of the samples, we evaluate A£ . In the A£ 
distribution for samples created without Bsi (^^2) 
highest observed value is 35.2 (41.8) in over 4.6 (5.6) 
million samples where the peak is located in range 0-50 
MeV/c^ (20 - 120 MeV/c^). Therefore, we conclude that 
the statistical significance of each of the signals exceeds 
five standard deviations. 

In summary, we report the first observation of the nar- 
row jq =3/2 states of the orbitally excited Bs mesons. 
The signals observed are attributed to the i?*2 — > B^K~ 
and Bsi B*^K~ decays. From the precise measure- 
ment of the Q values, we derive the masses of the two 
states and their mass difference, and the values are con- 
sistent with theoretical predictions [11 . 
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